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We present specific heat, cp, and ARPES data on single crystals of the stoichiometric super- 
conductor LiFeAs. A pronounced anomaly is found in cp at the superconducting transition. The 
electronic contribution can be described by two s-type energy gaps with magnitudes of approxi- 
mately Ai = 1.2 meV and A2 = 2.6 meV and a normal-state 7 coefficient of 10 mj/mol K 2 . All 
these values are in remarkable agreement with ARPES results. 



PACS numbers: 74.25.Bt, 74.25.Jb, 74.70.Xa 

A. Introduction 

The symmetry of the superconducting order parameter 
is one of the basic characteristics of the superconducting 
state. In this respect, the recently discovered pristine su- 
perconductor LiFeAs with a T c of approximately 18 K 
plays a decisive role in elucidating the pairing mecha- 
nism of the Cooper pairs and the nature of the supercon- 
ducting state in pnictide superconductors. In contrast 
to other pnictides, superconductivity in LiFeAs evolves 
without additional doping, and nesting between hole and 
electron pockets is very poor. The evolution of spin- 
density wave (SDW) type magnetic order, which is typi- 
cally present in the vicinity of the superconducting state 
in so-called '1111' and '122' pnictide superconductors, is 
not observed in LiFeAs. Consequently, mediation of su- 
perconductivity by antiferromagnetic spin fluctuations as 
suggested for other pnictides 5,6 is unlikely. Remarkably, 
there is evidence both from theory and experiments for 
almost ferromagnetic fluctuations which drive an insta- 
bility toward spin-triplet p-wave superconductivity. ' 

Experimental investigations on the structure and mag- 
nitude of the superconducting gaps in LiFeAs by means 
of bulk specific heat data as well as by angle resolved 
photoemission spectroscopy (ARPES) are of great in- 
terest. Previous specific heat data obtained on an as- 
sembly of tiny LiFeAs single crystals are in agreement 
with the presence of two isotropic gaps of 0.7 meV and 
2.5 meV, although the presence of nodes could not be 
ruled out. 11 Similarly, magnetization measurements on 
polycrystals suggest two s-type gaps of 0.6 meV and 3.3 
meV, but do not exclude gap nodes. 1 " Recent magnetiza- 
tion data of single crystals revealed two gaps of approx- 
imately 1.3 meV and 2.9 meV. 1 ' Measurements of the 
London penetration depth of single crystals are in line 
with nodeless superconductivity and two gaps of 1.7 meV 
and 2.9 meV. 12 ' 13 First ARPES measurements suggest 
the presence of two gaps of 1.5 meV and 2.5 meV as 
well 4 , with the larger one being in reasonable agreement 



with the results of the data analysis presented in 14. 

In this manuscript we study the superconducting en- 
ergy gaps of LiFeAs single crystals by two complemen- 
tary experimental methods, specific heat measurements 
and ARPES. The specific heat is sensitive to the bulk 
and gives direct access to the entropy of Cooper pair 
breaking, which is determined by the superconducting 
gap structure. In turn, ARPES allows for probing the 
momentum-resolved superconducting gap. From our re- 
sults we can exclude the possibility of d-wave pairing in 
LiFeAs. Instead, both methods are in line with the ex- 
istence of at least two s-type energy gaps for different 
Fermi surface sheets of LiFeAs with magnitudes of ap- 
proximately Ax = 1.2 meV and A 2 = 2.6 meV. 



B. Experimental methods 

Single crystals of LiFeAs have been prepared by a self- 
flux method described in detail in Ref. 15. Specific heat 
data were obtained by a relaxation technique in a PPMS 
(Quantum Design) on a sample with a mass of 2.4 mg. 
Measurements of the magnetic susceptibility in a field 
of 2 mT (MPMS-XL from Quantum Design) confirmed 
that this sample has a magnetic T c x of 16.9 K, similar 
to other crystals from the same batch. 1 ' After correction 
for demagnetization effects using an ellipsoid approxima- 
tion for the sample shape, the superconducting volume 
fraction was estimated to 0.91. The difference to 100 % 
is within the error of demagnetisation and may be due to 
deviations of the sample from an ellipsoid and/or due to 
a small non-superconducting phase. Photoemission ex- 
periments have been carried out using the synchrotron 
radiation from the BESSY storage ring. The end-station 
"1-cubed ARPES" is equipped with a 3 He cryostat which 
allows to collect angle-resolved spectra at temperatures 
below 1 K. All single crystals have been cleaved in UHV 
exposing the mirror-like surfaces. 
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FIG. 1: The temperature dependence of the specific heat of 
LiFeAs is shown as c p /T vs. T in zero magnetic field and in a 
field B || c of 9 T. The phonon and electron contributions c p h 
and Cei-n 9T have been determined from a fit of the 9 T data 
as explained in the text. In addition, the fitted total c p /T for 
T and 9 T are shown. 



C. Results 

The temperature dependence of the specific heat c p of 
a LiFeAs single crystal is shown in Fig. 1 as c p /T vs. T. 
In zero magnetic field a clear anomaly is observed around 
15 K which is attributed to the superconducting phase 
transition. By applying a magnetic field B = [j,qH \\ c 
of 9 T the anomaly is shifted to lower temperatures and 
reduced in height. 

In order to determine the specific heat related to the 
superconducting phase transition we need to estimate the 
phonon (c p h) and electron (c e i- n ) contributions to c p in 
the normal state. At low T, c c \_ n behaves linear in T, 
while Cph varies as c p h oc T 3 . However, for LiFeAs the on- 
set of superconductivity limits the fitting range towards 
low T. In order to improve the reliability at higher T, we 
used a second term of the harmonic-lattice approxima- 
tion, i.e. Coi-„ + Cph = 7T + /3 3 T 3 + /3 5 T 5 . The results of 
a fit of the 9 T data between 13 K and 20 K are shown as 
lines in Fig. 1. From the fitting parameters we calculated 
a Debye temperature of 310 K and a Sommerfeld coeffi- 
cient of 7qt = 8.5 mJ/mol K 2 . Below 13 K the 9T data 
deviate slightly from the fit, although the superconduct- 
ing transition is observed only around 8K. This may be 
due to a tiny part of the sample with different orienta- 
tion B || ab. For this field direction the superconducting 
transition in 9 T is known to take place at 13K. 16 The 
data above 13 K are well described by the fit. 

The zero field specific heat is shifted by a constant 
with respect to the 9T curve above T c , which is at- 
tributed to a field-dependent electronic specific heat. 
From our data we obtain a zero-field Sommerfeld coef- 
ficient 7qt = 7n = lO.OmJ/mol K. This value is about 
half the one determined in previous specific heat studies 



on polycrystalline LiFeAs and on an assembly of small 
single crystals/' The differing 7 coefficient may arise from 
differences in the sample quality. Resistivity measure- 
ments on a crystal from the same batch as the one investi- 
gated here revealed a very low residual resistivity of only 
15.2 //ficm and a large residual resistivity ratio RRR = 
P300k/Pok of 38. 18 This confirms the high quality of our 
samples. By contrast, the polycrystalline LiFeAs investi- 
gated in Ref. 17 exhibits a much larger residual resistivity 
of approximately 2.5 mf2cm with RRR w 10. This may 
be an indication for the presence of impurities, which can 
give rise to additional contributions to the specific heat. 

Fig. 2 shows the temperature dependence of the elec- 
tronic contribution to the specific heat c e \ in zero field 
determined by subtracting c p h- From an entropy-con- 
serving construction we find a superconducting transition 
temperature of T c = 14.7 K, which is somewhat smaller 
than the magnetic T c x of the sample. This difference 
between the thermodynamic and magnetic T c has been 
also found in previous studies of LiFeAs. 1,1 Of partic- 
ular interest in this context is a recent investigation of 
single-crystalline LiFeAs: 11 ' The onset value of T c deter- 
mined by dc susceptibility in 1 mT was found to be about 
17 K. However, the extrapolation of the field dependence 
T C (B) determined by magnetic torque measurements on 
a small piece from the same crystal yielded a zero-field T c 
of 15.5 K. Likewise, ac-susceptibility data for B > 0.5 T 
on crystals from the same batch extrapolate to a lower 
T c than the values measured in very small fields. There- 
fore, the difference between the low-field T c x and the bulk 
T c determined, e.g., by specific heat appears inherent to 
LiFeAs. 

The jump height of c c \/T at T c amounts to 1.24 7„, 
which is slightly lower than the BCS value of 1.43 7„. 
In addition, the almost linear temperature dependence 
of Cei/T indicates, that the specific heat data cannot be 
described by a single BCS gap. In order to illustrate this 
we show a theoretical BCS curve with A = 1.764 k-B,T c = 
2.23 meV in Fig. 2, where we account for a small frac- 
tion 7 r es/7n = 0.09 of normal electrons as justified be- 
low. Systematic deviations of the single-gap fit from the 
data are observed in the whole temperature range below 
T c . Since a single gap cannot describe the data, we ap- 
plied a phenomenological two-gap model developed for 
the specific heat of MgB 2 . j9 For this purpose we cal- 
culated theoretical curves c c \/T vs. T for a large range 
of the free parameters Ai, A2, 7 rC s, and the weight of 
the gaps wa2/ai- We then calculated the differences di 
between these curves and the data point by point and 
used the sum of the squares £d 2 as criteria for the model 
quality. The best description of the data is obtained for 
Ai = 1.44 meV, A 2 = 2.74 meV, 7 rcs /7„ = 0.13, and 
w A2/ai — 1-18. This result, however, is not realistic: 
The measured data exclude values larger than 0.11 for 
7ros/7n- This is seen from the inset of Fig. 2, which shows 
the low-temperature ratio c c \/^ n T. It reaches a value of 
0.11 around 2 K. In view of the systematic decrease of 
the data with decreasing temperature it is likely, that the 
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FIG. 2: The main plot shows the electronic contribution to the 
specific heat of LiFeAs in zero field as c c i/T vs. T. The dotted 
line corresponds to a single-gap BCS curve taking into account 
a residual nonsuperconducting contribution with 7 rcs /7n = 
0.09. The red curve is a two-gap model of the data assuming 
likewise 7 rcs /7 n = 0.09. The inset shows the low-T data as 
CeiHnT vs. T on a larger scale. In addition, an estimate for 
the low-T Cei expected in case of a d-wave order parameter is 
plotted. 



ratio 7r CS /7n is even lower. Therefore, we take the normal 
state contribution determined from the magnetic suscep- 
tibility of the same sample to estimate 7 rcs /7n = 0.09. 
By doing so, the best description of the data is obtained 
for Ai = 1.24 meV, A 2 = 2.62 meV, and w A 2/ai = 1-53. 
The corresponding calculated specific heat is shown as a 
red line in Fig. 2. It is in a very good agreement with the 
data in the whole temperature range. Still, a close look 
to the low-T part (cf. inset of Fig. 2) reveals systematic 
deviations from the data, which suggests an even lower 

7res- 

Although the data are best reproduced by the parame- 
ters given above, they can be described with similar accu- 
racy for a considerable range of gaps. For an estimate we 
consider all curves for which the deviation from the data 
Edf is at most 2 times the value for the best curve with 
7rcs/7n = 0.09. As an additional constraint we assume 
7rcs/7n < 0.11 in agreement with the data. Thus, we ob- 
tain Ai = (0.93 - 1.67) meV, A 2 = (2.40 - 3.24) meV, 
7r CS /7n = 0.04 - 0.11, and w A2 /ai = 0.46 - 3.45. 

The overall shape of the superconducting anomaly pre- 
sented here is similar to the one obtained recently on a 
cluster of tiny single crystals of LiFeAs.'' However, the 
magnitudes of the gaps are somewhat larger, in particu- 
lar for the smaller gap of 1.24 meV compared to 0.7 meV 
obtained in Ref. 9. This may be due to the limited reso- 
lution of the data in Ref. 9, which leaves a considerable 
uncertainty for 7 rcs . In addition, the magnitude of the 
smaller gap Ai itself is very sensitive to the specific heat 
at low T. The high resolution of our data down to 2 K 
allows for a very reliable estimate of both 7 ros and Ai . In 



FIG. 3: (a), (b) Photoemission intensity in the high sym- 
metry cuts together with the fitted tight-binding dispersions, 
ei, e2 are the two bands forming the electron pockets at the 
X point of the BZ. hi, I12 are the outer and the inner hole 
bands centered at the Y point, (c) Band contribution to the 
DOS derived from the fitted quasiparticle dispersions, (d) 
Temperature dependent Sommerfeld coefficient determined as 
7 (T) = dc cl /dT. 



Ref. 9 even a d-wave order parameter could not be ruled 
out. In this case one should find a linear-in-T behavior 
of c c \/T for T <C T c . Since our measurement range is 
limited to T > 0.14 T c , we cannot directly exclude the 
presence of line-nodes of the gap function, either. How- 
ever, the quality of our data and the low value at 2 K ren- 
der the presence of line nodes of the gap function very 
unlikely. This is demonstrated by the expected low-T 
behavior for a d-wave order parameter estimated as de- 
scribed in Ref. 9 and assuming 7 rcs = 0. It cannot be 
brought into agreement with the measured data. Taking 
into account in addition recent ARPES results suggest- 
ing a nearly isotropic gap for each Fermi surface sheet, 4 a 
d-wave order parameter is excluded from consideration. 

Owing to its ability to resolve both momentum and en- 
ergy of the electronic states, ARPES can provide a com- 
plete picture of the electronic band dispersion and the 
momentum-dependent superconducting gap. Therefore, 
it is interesting to compare the thermodynamic proper- 
ties measured directly with those calculated from photoe- 
mission data. In particular, one can quite easily extract 
the value of the Sommerfeld coefficient 7 n . This param- 
eter determines the heat capacity in the normal state, 
and, along with the superconducting gap values and the 
assumption of BCS pairing, the thermodynamic proper- 
ties in the superconducting state. For this purpose the 
photoemission intensity of the LiFeAs has been mapped 
in more than one Brillouin zone (BZ) and fitted with 
the standard tight-binding formula. 20,21 To demonstrate 
the agreement between the ARPES raw data and the ob- 
tained tight-binding fits, in Fig. 3a and b we show two 
high symmetry cuts along T-X-T and T-M-r directions 
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with the fits to the renormalized quasiparticle dispersions 
superimposed over the ARPES data. As one may see 
from the derived density of states (DOS) (Fig. 3c), the 
major contribution to the heat capacity must be due to 
the outer hole band. Another interesting observation is 
a quite pronounced variation in the DOS at the Fermi 
level, which may result in deviations from the linear tem- 
perature dependence of the electronic heat capacity. To 
check to which extent this applies to the current case, 
instead of using a standard text-book expression for the 
electronic heat capacity c c \ oc V{E-p)k-QT ', we have made 
an estimate based on a more general expression, which 
for a single quasiparticle band with dispersion £ ki reads 
as c M = 2N A -^ r (f(£ ki ,T)£ ki ) B z- Here (...) B z denotes an 
average over the Brillouin zone and N\ is the Avogadro 
constant. Indeed, the value of dc c i/dT, that in case of 
strict temperature linearity defines the Sommerfeld coef- 
ficient 7„, varies from 13 to about 17mJ/mol K 2 , which 
is in a relatively good agreement with the direct measure- 
ments resulting in 7„ ss 10 mJ/molK 2 . This variation in 
the dc e i / dT is also likely to account for some variation in 
the 7 n extracted from the thermodynamic measurements 
by different authors. 

It is remarkable that the thermodynamic gap values 
determined in the present study are in excellent agree- 
ment with the ARPES leading-edge gaps reported in 
Ref. 4 (1.5 meV and 2.5 meV). However, the absolute 
values of the actual gaps, which can be derived from the 
ARPES data after a more rigorous analysis, are usually 
slightly higher than the leading-edge ones. The resulting 
discrepancy (of the order of 0.5 meV) between the abso- 
lute gap values derived from ARPES and in the present 
study could probably be explained by the difference in 
thermodynamic and magnetic T c mentioned before. In 
any case, the ratio of the ARPES leading-edge gaps re- 
produces the ratio of the actual gap values quite accu- 
rately, and this is in close correspondence with the ones 



discussed here. Moreover, a more detailed investigation 
of the superconducting gaps in LiFeAs 22 indicates that 
the value of the gap supported by the band h2 is in- 
deed comparable to that of the large hole Fermi surface, 
made by the band hi (see Fig. 3a). In the light of the 
results presented in Fig. 3c, it is important to establish 
the fact that the smaller gap corresponds to the hole-like 
Fermi surfaces centered around the Gamma point, while 
the larger one corresponds to the electron-like Fermi sur- 
faces localized around the corners of the Brillouin zone. 
This knowledge can help to identify the symmetry of the 
order parameter in iron pnictides in more details. 



D. Summary 

In summary, both thermodynamic and spectroscopic 
experiments on LiFeAs render a nodal gap very unlikely, 
and equivocally speak in favor of a strong variation of the 
gap magnitude between different electronic bands, from 
about 1.2 meV to 2.6 meV. The general agreement of such 
complementary probes within band picture emphasizes 
the robustness of the conclusions drawn. The multigap 
behavior of LiFeAs established above, is in line with two 
gaps found in many other iron arsenides. 
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